Equations in dual variables for Whittaker functions. 
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Abstract 

It is known that the Whittaker functions w(q, A) associated to the group SL(N) are eigen- 
functions of the Hamiltonians of the open Toda chain, hence satisfy a set of differential equations 
in the Toda variables qi. Using the expression of the qt for the closed Toda chain in terms of 
Sklyanin variables A^, and the known relations between the open and the closed Toda chains, 
we show that Whittaker functions also satisfy a set of new difference equations in A^. 
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1 Introduction. 



The method of separation of variables is acquiring a central place in the domain of integrable 
systems. In the classical case, the separated variables are simply the poles of the eigenvectors of 
the Lax matrix. Together with the spectral curve, they are the necessary data to reconstruct the 
eigenvectors from their analyticity properties, and therefore the Lax matrix itself. The solution 
of the integrable system is given the fact that the image of the divisor of poles by the Abel map 
evolves linearly on the Jacobian under an integrable flow. 

In the quantum case, these separated variables where first defined and used by Sklyanin. In partic- 
ular, he showed that the original iV-body Shroedinger equation separates into N one dimensional 
equations, called Baxter equations. The separated variables also appear in the iV-soliton form factor 
formulae in sine-Gordon theory. It was shown that Smirnov formulae have a simple interpretation 
in terms of the separated variables of the N solitons [7j. 

More recently, it was realized that the separated variables provided a convenient approach to 
the problem of quantizing a classical integrable system. It is interesting that this idea emerged 
independently both among mathematicians JJ| and physicists |12j . 

In [Tl], we took this point of view to solve the quantum inverse problem for the closed Toda 
chain. In this note, we remark that this result immediately suggests that the Whittaker functions 
associated to Sl(N), which are the eigenfunctions of the open Toda chain, see eas. (fTT|) below, should 
also satisfy a set of equations in the dual (momentum) variables, eqs. (|18|) below. We prove these 
relations by using the Mellin-Barnes integral representation for the Whittaker functions found in 

EH- 



2 The Closed Toda Chain. 

The closed Toda chain is defined by the Hamiltonian 

N+l 
i=l 

where we assume that qN+2 = 1i, an d canonical Poisson brackets 

Qj} = 0, {Pi, qj} = <%, {pi,Pj} = 



As it is well known, one can associate to this system a 2 x 2 Lax matrix as follows. Consider the 
matrices 

A + pj -e q i 
e~* 



Tj(X) 
and construct 



T(X)=T 1 (X)---T 2 (X)T N+1 (X) (2) 
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We can write 

T(A) = (c(A) d{\)) ; A(WA)-J3(A)C(A) = 1 (3) 

A(X) is a polynomial of degree N + 1, -D(A) is of degree iV — 1, and B(X), C(A) are of degree n. 
The spectral curve is denned as usual 

det(r(A) - fx) = = // + - t(A) = (4) 

where 

* 1 
t(A) = ^(A)+D(A) = A Ar+1 + ^A^ J , ffjv = -P, H N ^ = -P 2 -H 

j=o 

where P = ^Pi, and H is given by eq.Q. We will consider the system reduced by the translational 
symmetry. We fix the total momentum P = 0. The symplectic form uj = J^j^ 1 dpi A dqi becomes 
^reduced = Y,i=\dpi A d(qi - qN+i) so that the canonical coordinates of the reduced system can 
be taken as (pi,qi — ?jv+i) 3 £ = 1---N. We choose the gauge condition qN+i = 0. It must be 
emphasized that this reduced system, which does not contain the degree of freedom pn+i,Qn+i 
anymore, is not the open Toda chain. The N quantities Hj are conserved and Poisson commute. 

The separated variables are the poles of the eigenvectors of the Lax matrix. For a 2 x 2 matrix of 
the form eq.© the eigenvector is simple 

to -,)•=* • = * = -^ 

The poles of ip2 a t finite distance are above the zeroes Aj of -B(A) = which is a polynomial of 
degree n. The two points above Aj are fif = A(Xi), fi~ = D(Xi) so that ^2 has a pole only at the 
second point. The points of the dynamical divisor are therefore (Aj,-D(Aj)), -B(Aj) = 0. 

The Inverse Problem consists in reconstructing everything in terms of the (Aj,/ij). Its solution for 
the quantum theory was given in |14j . 

Quantum commutation relations are defined directly on the separated variables. 

[Xk,Xk>] = 0, HkXk> = (Afc/ + ih5kk')^k, [fJ>k,Hk>] = 

To reconstruct the Hamiltonians themselves is simple. The n points Aj, /i~ belonging to the spectral 
curve, we have the equations (we drop the superscript — in fj,~): 

N 

w+^ 1 -Af +1 -E A ^ = ° 

3=1 

This is a linear system of N equations for the iV quantities Hj. Its solution makes sense even 
quantum mechanically and, quite generally, produces commuting Hamiltonians [111 I12j . 

To reconstruct the original operators Pi,qi, we need the operators X^\Y^ below. We call [k] a 

subset of cardinality k of (1, 2, • • • , N): [k] = (£1, £2, • • • , £&)• We write Yl[k] for the sum over all 
such subsets. Let 

S lk]( X ) = II II (A -A V ^l =ftlft2 '"^ ( 5 ) 
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and define 



X(k) = 52S [k] (\) m (6) 

[k] 

Y(k) = £%](*)(£*] mm ( ? ) 

[k] um J 



The quantum Toda operators are given by: 



" l J / ST ^ ' 



The canonical commutation relations 

[ e *, e %] = J [e*,pj] = -ihS ije qi , \pi,Pj}=0 (9) 
are a simple consequences of the following quadratic algebra: 
[XW,X®] = 

[X {k \Y^\ = ih(k - l)6(k - l)X^ x® 

[y (k) 5 y (0] = ^ _ /) r fl ( fc _ j)y(*)jf(0 + 0(/ _ ijyWxW" 

Note that this implies that there is no ordering ambiguity in the expressions eq.(|5]>. 

The operators are self adjoint with respect to the non trivial scalar product 

/oo 
d\m{\) f*(X)g(X) 
-oo 

where 



m 



3 Whittaker vectors and functions. 

There is an important connection between the closed and open Toda chains that we now recall. 

The matrix element B(X) in eq.© can be written as 

N 



k=0 

where hk{p,q) are the Hamiltonians of the open Toda chain obtained by removing particle N + 1 
from the closed chain 2 . 

The eigenfunctions of the Hamiltonians hk(p,q) satisfy 

hk{p, q)w(X, q) = a k (X)w(X, q) (11) 



where cr^(A) are the symmetric functions of the Aj. It is a beautiful result by Kostant and Semenov- 
Tian-Shansky that w(X,q) is the Whittaker function of sl(N). Thus the function w(X,q) is the 
kernel of the Fourier transform going from the q to the A (momentum) variables. They satisfy the 
completness relation 

r dX m(X)w*(q, X)w(q', A) = 6(q - q') (12) 



where m(X) is given by ea. (|l(J|) . and the orthogonality relations [5] 

J dq w*(q, X)w(q, A') = m _1 (A)£(A - A') (13) 

To recall how such functions arise, consider the sl(N) Lie algebra 

[Hj, E± ai ] = ±aijE± ai , [Ecu, E~ aj ] = dijHi, i,j = l,---,N-l 
where on are the simple root vectors and a,j is the Cartan matrix. The Weyl vector is 

The quadratic Casimir operator is 

a i,j 

where the ^2 a runs over all positive and negative roots. 
Whittaker vectors are such that 

E ai \wx) = Vai\ w x)i a i simple (14) 

Of course, if a is a non-simple root, then E a \w\) = 0. Similarly, dual Whittaker vectors satisfy 

(w\\E-a. = ^{wxl, ai simple (15) 

We may furthermore assume that \w\) belongs to some irreducible representation with weight A. 
In that case, we have 

C 2 \w\) = c 2 (X)\w\) 

Define the Whittaker function 

w(X,q) = e p{q) (w x \e q \w x ) 
where q = ■ qia~^Hj belongs to the Cartan subalgebra. Then, we have 

c 2 (X)w(X,q) = e^(w x \e q C 2 \w x ) 

\ a>0 i,j J 

Expanding H a = ^j( a )« a ij 1 ^i' wnere a ij is the Cartan matrix, we get 

C2 (X)w(X,q)=L,^- + 2 J2 ^- a(9) U(A,g) 

\ 1 3 a simple / 
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The differential operator in the right hand side is just the Hamiltonian of the open Toda chain. 
Hence w(\, q) satisfy to the open Toda chain Schroedinger equation! Clearly, the same analysis can 
be done for higher order Casimirs giving rise to the higher order Hamiltonians of the open chain. 



There exists several integral representations for the Whittaker functions. 



The "Gauss" representation jOj. 



w G 



. N-l 
(X,q) = / dZij JJ A 
i=l 



The integral is over upper-triangular matrices N x N matrices, z, with 1 on the diagonal. S 
is the antidiagonal matrix S{j = $N+i-i.j- Aj(M) is the determinant of the i x i submatrix 
of M consisting in the first i rows and columns. Aj_i ; j(M) is defined as the determinant of 
the (i — 1) X (i — 1) submatrix of M with columns i — 1 and i interchanged. 

The "Iwasawa" representation jHj. 



wi(X,q) = e 



Kg) 



r N ~ l i i 
J i=i 



The definitions of z and Aj(M) are the same as above. 
The "Mellin-Barnes" representation [TU] . 

The Weyl invariant Whittaker function has a representation in terms of multiple Mellin- 
Barnes integrals. Let ■jjk be a lower triangular N x N matrix. 

o ... \ 

722 



/ 7n 

721 



7Af-l,l TJV-1,2 
V 7JV1 1N2 



TfN-lytt-1 

JNN / 



We identify 
Then 



ry Nj = \ j) j = l---N 



N-l 



v-i fVr n tt"+i f, lni ~ 1 r n+Uk r fy-:-- „ u. 
lij=iiik=i n lh L \ m 



(16) 



n 



yin ir/ 7nj-7n,t M2 
n=l ^ llj<k=l \ L \ ih )\ 

The integration contour is as follows. 

Im 7n > max(Im 712, Im 722) 
Im 721, Im 722 > max(Im 731, Im 732, Im 733) 



Im 7jv-i,i, • • • ,Im 7iv-i,7V-l > max(Im 7^1, • • • , Im 7^) 
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Note that the poles of the integrand are located at 

j nj = 7n+i,fc - iHs, 
so we can move the contour j n j upward safely. 



4 Diagonalization of X^ k \ 

Since the operators are all commuting, we can diagonalize them simultaneously. The X^ are 
self-adjoint with respect to the measure m{\) eq.QUJ), their eigenfunctions satisfy the same com- 
pleteness and orthogonality conditions as the Whittaker functions diagonalizing the Hamiltonians 
of the open Toda chain eqs. (|12J13|) . 

The Whittaker functions are the kernel of the Fourier transform going from the variables qi to 
the variables \. Similarly, the eigenfunctions of the opertors X^ are the kernel of the Fourier 
transform going from the variables Aj to the variables q{. So, it is natural to expect 

X( fc )(A ;/ u) w*(\,q) = (-l) k (N-i) e ELiiN + i-,. w *(\ : q) (17) 

Taking the complex conjugate of this equation, we get 

X«*(A, m) w(X, q) = (-l)^" 1 )^! <iN+i-i W ( A) q ) (i 8 ) 

where X^* is the complex conjug ate of i.e. it is given by eq.© with [ii replaced by \x\ the 

shift operator \ — > Aj — ih. Note that eas. ()ilj) are differential equations in qi, while eas. (|TH|) are 
difference equations in Aj. 

In the remaining of this section, we prove that the function w(\,q) defined by the Mellin-Barnes 
integral representation eq. (fTK|) . which is known to satisfy eqs. (fTTj) . also satisfy eqs. ()18|) . 

Before treating the general case, it is instructive to do the calculation for N = 2 and N = 3 first. 



4.1 The case N = 2. 

For N = 2, the function u>(A, q) reads: 



The operators X^* are 



xW* = -±—^1-^), x^'=^* 2 

Ai — A2 
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We have 
X {ir w(X,q) = 

(711 - Ai) (Til - ^2) 

; "" 1 V — 77^ — ^ 

k=l 



ih ih 



2 , . 

~<ii- x k I 711 - X k 



^ f d 7 ll e -mfe7n+ g2 (A 1 - 7 n + A 2) ] rj (a) ^ r 

12 J fc=1 V z / 

r 2 -a 

-e 92 / d7n e -m[9iTu+?2(A 1 -7ii+A 2 )] Y[(ih) m ^T . 

fe=i ^ 1 

Hence, we have proved 

X^*w(X,q) = -e q2 w(X,q) 

Similarly 

X(»w(\,q) = J d 7ll e-^i+^m+x^m)) -Q ( .^m^ r ^ 7u - A fc + ^ 



where 7^ = 711 + ih. So, the integral on the right hand side is the same as the initial one but on 
a contour shifted upward by +ih. But as we already noticed, this does not change the value of the 
integral. Hence, we have proved 

X( 2 >w(\,q) = e qi+q2 w(X,q) 



4.2 The case N = 3. 

For N = 3, we have: 

w(X,q) = / ^ e -i7i [91711+92(721 -7n+722)+iJ3(Ai+A2+A 3 ~72i -722)] x /-j_g^ 



X 



nip / 7lj-7ifc \ 12 TT Ip / 72j~72fc \ 12 

i<fe I 1 I m I I iij<k I 1 I in 



The operators X^* read 



A12A13 A21A23 A31A32 

X = 1 \ — ^V" 2 + \—^—VlV3 + 1-^-/^3 

A13A23 A12A32 A21A31 

Consider first X^*w(X,q). The exponential factor in the first line of ea, (|19j) produces the factor 
e qs . In the integrand, using T(x + 1) = xT(x), we get the factor 

(ih) 2 (721 - Ai) (722 ~ Ai) (ih) 2 (721 - A 2 ) (722 ~ A 2 ) (^) 2 (721 ~ A 3 ) (722 - A 3 ) 
A12A13 ih ih A21A23 ih ih A31A32 ih ih 
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which is equal to 1. Hence, we have proved 

X^*w(X,q) = e q3 w{X,q) 



Next, we consider X^*w(X, q). The exponential factor produces the factor e 2<?3 . Using the T 
function relation, we get in the integrand 



A13A23 
1 

A12A32 
1 

A21A31 



(721 


- Ai)(722 


- Ai)(72i 


- A 2 )(722 


-A 2 ) + 


(721 


- Ai)(722 


- A1X721 


- A3) (722 


-A3) + 


(721 


- A 2 )(722 


- A 2 )(72i 


- A3) (722 


-A3) 



which is equal to 

(721 - A1K721 - A2X721 - A 3 ) - (722 - Ai)(722 - A2X722 - A 3 ) 

(721 - 722) 

So, we get a sum of two identical integrals but with 721 and 722 interchanged. Let us treat the first 
one. Using the T functions, we reconstruct 



r|77i-o^ix 



n fc r(^ + i)r(^ 

ih ) r ( 721 ^ 722 + 1) r ( 722 ^ 721 ) 

We now change variables 721 = j' 2 i — ih. The exponential produces the factor e q2 ~ qs . The integrand 
becomes 

TT ]-> ( 721 -Afc \ p / 722 -A fc 
I . , 711 ~ 721 \ p / 711 ~ 722 \ V ih ) \ ) 711 - 721 



ih J \ ih J r( 721 ~ 722 )r( 722 ~ 721 ) 722-721 

Note that there is no pole at 721 = 711 nor at 721 = 722; so that we can move back the 712 contour 
to its original position. To this, we have to add the same expression with 721 and 722 interchanged. 
We reproduce w(X, q) because 

7n ~ 721 7n — 722 _ j 
722 — 721 721 — 722 

Thus we have proved 

X^*w(X,q) = e q2+qs w(X,q) 

Finally, we consider X^*w(X, q). The exponential factor produces the factor e 3q3 . Then we change 
variables 72^ +ih = 7^, 711 + ih = 7^, which amounts to shifting the contours by ih upward. This 
produces a factor e qi+q2 ~~ 2q3 . So we have proved that 

X^*w(X,q) =e qi+q2+q3 w{X,q) 
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4.3 General case. 



We have 

X( k >w(\,q)=Y,S [k] tf k] w(\,q) 

[k] 

The exponential in the first line of ea.(|16|) produces a factor e kqN . In the integrand we get the 
factor 

N-l 

Ik] je[k] »=l 

By eq. 1)20(1 in the Appendix, this is equal to 

N 

(_!)(*-*)* j2 ^-1) n n^-M-T^-) 

[fe-i] ie[fc-i]i=i 

Let us keep track of 7jv-i,ij i £ [& — 1]. The factors n«e[fc-i] rij=i(7iV-l,i — 7N,j) are absorbed into 
the T-functions in the numerator of the integrand, while the factor Sjfc-i] (72V— l) is absorbed in the 
T-functions in the denominator of the integrand, to produce 



n v r 1 7Jv-i,i-7Arj , 1 



ie[fe-i] llj0[fc-i] 1 ^ il h J- J 1 ^ m 



Next, we change variables jN—l,i + ih = 7yv-i j, i G The exponential factor in the first line 

of eg. 1)16)1 is symmetrical in all the 7jv— l,i an d yields 

e (fc-l)(gjv-i-giv) 

Then, the T-functions with shifted arguments are 

tN-2 -p ( lN-2,j-lN-l 



T-rA t-2 y I lN-2,j-lN-\,i , ^ 



n W 7jv-i,j-7jv-i,j 1 1 



ie[fc-l] llj£[fc-i] 1 ^ W 
which produce a factor 

N-2 

(-i)( iv - fc )s [A; _ 1] (7^i) n n^-2 J -7iv- M ) 

ie[k-i] j=l 

in the integrand. So, we are back to the same problem, but at level k — 1,N — 1. We eventually 
reach the level 1, N — k + 1, where we use the identity 

TV-fc 

£ n w . w - TO . tHJ n - - (-1)'-" 

Putting everything together, we arrive at 

X (k >w(X,q) = (-l) k ^ N -^e qN+qN - 1+ - +qN ~ k +'w(X,q) 
which is exactly eq. 1)18(1. 
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5 Conclusion. 



To conclude, we would like to mention that the operators are limiting cases of the Ruijsenaars- 
Macdonald operators: 



tq Xi - t _1 g A J 



^ = e n -h^ 1 ^ 



[k] »e[fc] 

jew 



So, our result is probably a limiting case of the results obtained in [P3] . 

Acknowledgements: I thank F. Smirnov for discussions. This work was partially supported by 
the European Commission TMR program HPRN-CT-2002-00325 (EUCLID). 



6 Appendix. 

Proposition 1 One has the identity 

N-l N 

£ 5 [*](™)II n^-W-Wi) = (- 1 ) (iV " fc) *E 5 [*-i]^-i) II E[(W-W-7W) (20) 

[k] je[k] i=l [k-l] ie[k-l]j=l 



Proof . First of all, the left-hand side is a polynomial in 77V. It has potential poles at jNi = iNj- 
Let us suppose i = 1, j = 2. One has to assume that 1 £ [fc],2 ^ [k] or vice versa, otherwise there 
is no pole. In the above sum, we consider the two terms 

[k} = l + [k% [k} = 2 + [k'}, l,2$[tf] 

Denote by [I'] the complementary subset of [k'] in 3, 4, ■ ■ • , N. The two terms are, respectively 

1 T — r 1 T — r 1 t — r 1 ^ ^ 



n — z — : n — ~z — n — ~ — : n (tw-u - im) n (tiv-u - im) 



INI ~ 1N2 INI - INj 77V, " 77V2 7iVi " 7iVj ^ 



1 rr 1 rr 1 rr J 



JV-1 



n n n n (w-u - 7^2) n - 7^0 

J6[I'] 

It follows that the sum of the two residues cancel. From the behaviour at oo, we see that both 
sides are polynomials of degree k — 1. To show that they are identical, we compare the values at 
the N — 1 points 7^1 = 7iv-i,v It is enough to consider 

7iVl = 7iV-l,i 

In the left hand side, only the sets [k] such that 1 [k] contribute. Hence we get 

1 1 N ~ 1 

e n lm - lm n 7/v? ._ 7/v , n - w n - w 

i2[fc]ie[fc] lNl lN1 *€[*] 7m 77Vj je[fe] i=2 
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When evaluated at 7jvi = 7iv-i,i, the first and third terms cancel and we are left with 

N-l 



E II zi:IIII(^-7«) 



mk] .JfW lNi 7Nj ^=2 je[k] 



In the right hand side, only the sets [k — 1] such that 1 ^ [k — 1] contribute. Hence, we get 

(-!)("-*>*£ n — i — n — - — n ^-m-t^iW-m-^) 

When evaluated at 77V i = jn— 1,1, the first and third terms cancel and we are left with 

( _i)C"-*)* y: n — - — n iW-^-^) 

The two things are identical if our identity holds at level N — 1. The lowest level is N = k. There, 
we have 

S[N](1n) = 1, «S'[jv_i](7jv-i) = 1 

and the identity reduces to 

N N-1 N-1 N 

n n (7iv-i,i - iNj) = yi n(jN-i,i - iNj) 

j=l i=l i=l j=l 

which is obviously true. ■ 
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